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One of the most important tasks in high energy physics is the precise 
determination of the strong coupling a s (Af|). Not only does measurement of 
a s (Af|) in different hard processes and at different hard scales Q provide a 
fundamental test of the theory of strong interactions, Quantum Chromodynamics 
(QCD), but it also allows constraints on extensions to the Standard Model of 
elementary particles [f]. The large set of a s (Af|) measurements is consistent with 
a central value of about 0.117 with an uncertainty of ±0.005 [2]. However, nearly 
all measurements are limited by theoretical uncertainties that derive from lack of 
knowledge of higher-order perturbative QCD contributions, or of non-perturbative 
effects, or both. It is hence vital to reduce the size of the limiting theoretical 
uncertainties which may, or may not, be concealing new physics. 

Here hadronic event shape observables in e + e - annihilation are considered. For 
an infra-red- and collinear-safe observable X: 

-%{X, n) = — s (n)A(X) + TilV) B(X, //) + al 3 (/i) C(X, n) + 0(af(fi)) (1) 

where = a s /2TT and // is the renormalisation scale. To date only the coefficients 
A(X) and B(X,fi) have been calculated [3,4]. The 15 hadronic event shape 
observables used in the recent a s (Af|) determination by the SLD Collaboration [5] 
were employed. Distributions of these observables were measured [5] using a sample 
of about 50,000 hadronic Z° decay events. The data were corrected for detector bias 
effects such as acceptance, resolution, and inefficiency, as well as for the effects of 
initial-state radiation and hadronisation, to yield 'parton-level' distributions, which 
can be compared directly with the QCD calculations. The EVENT program [6] was 
used to calculate the coefficients A and B in Eq. (1). 
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First, the 0(a 2 ) calculation, using the physical scale \i = Mz, was fitted [7] to the 
measured parton-level distributions by minimising % 2 w.r.t. variation of A^r. Each 
resulting A-^r value was translated into a s (M z ) and is shown, with experimental 
errors [5], in Fig. 1(a). There is considerable scatter among the 15 a s (M z ) values. 
Similar results have been observed previously [8]. Since the same data sample was 
used to measure each observable, and since the observables are highly correlated, 
this scatter is very significant and can be interpreted as arising from uncalculated 
higher-order perturbative QCD contributions, which a priori may be of different 
sign and magnitude for the different observables. The average a s (M z ) value and 
corresponding r.m.s. deviation are listed in Table 1. 

This procedure was repeated [7] using the experimentally-optimised-scale 
approach [9] in which a simultaneous fit of Aj^g and // to each observable was 
performed. Each resulting A-^-gr value was translated to a s (M z ) and is shown 
in Fig. 1(a). For -D-f no minimum in \ 2 w.r.t. variation of // in the range 
10" 4 < IJ 2 /M 2 Z < 10 2 could be found [7]. A gain, there is large scatter among the 
14 a s (M z ) values. For most observables the experimentally-optimised scale yields a 
lower value of a s (M z ) than the physical scale because the optimised scale is typically 
smaller than Mz, requiring a smaller value of Aj^g in order to fit the data [10]. The 
average a s (M z ) value and r.m.s. deviation are listed in Table 1. The r.m.s. deviation 
is comparable with that resulting from the choice of the physical scale, implying 
that use of the experimentally-optimised scale does not serve to reduce uncalculated 
higher-order effects. 

From the //-dependence a 'renormalisation scale uncertainty' on a s (M z ) can 
be defined [5] for each observable; these are shown as bands in Figs. 1,2. Within 
such uncertainties the a s (M z ) values determined from the different observables 
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using either the physical or experimentally-optimised scales are consistent, but this 
arbitrary procedure leads to a large uncertainty of ±0.0106 on the average value of 
a s (M 2 z ) [5]. 

The best resolution of this situation would be to calculate the observables 
to higher order in perturbation theory, a difficult and unattractive task that has 
not yet been achieved. In the absence of O(af) QCD calculations it has been 
suggested that the O(o^) calculations can be 'optimised' by choosing a specific value 
of the renormalisation scale. Since the all-orders result would be independent of 
renormalisation scale, Stevenson suggests that // be chosen according to the 'Principle 
of Minimal Sensitivity' (PMS) [11], so that da(X, fi)/dfi = 0. Grunberg suggests 
that fi be chosen to give the 'fastest apparent convergence' (FAC) of the series [12], 
so that the second-order term in Eq. (1) vanishes. Brodsky, Lepage and Mackenzie 
advocate that // be chosen to remove the iVy-dependence of the second-order term in 
Eq. (1), effectively incorporating quark and gluon vacuum polarisation contributions 
into the definition of the strong coupling [13]. 

For each observable the PMS, FAC and BLM optimised scales were calculated 
[7] and used in turn in a fit of the O(o^) calculation to each measured distribution 
to determine A^t and hence a s (Af|). The results are shown in Fig. 1(b); in 
the case of the oblateness O an acceptable fit with the BLM scale could not be 
obtained. For each observable the PMS- and FAC-derived a s (Af|) values are 
very similar, whereas, in some cases, the BLM-derived a s (Af|) value differs from 
them. This behaviour follows from the correlation between the scale value and the 
corresponding A-^r required to fit the data [10]. For a given observable the PMS- and 
FAC-derived a s (Af|) values are often, though not always, close to that determined 
using the experimentally-optimised scale. Furthermore, for most observables the 
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PMS-, FAC- and BLM-derived a s (M^) values all lie within the range encompassed 
by the renormalisation scale uncertainty defined in Ref. [5], though for p } Byy, B^-, 
D%°, Df and (B T ), the BLM- (PMS/FAC-) derived values lie below this range. 

For any of the PMS, FAC or BLM scale choices there is considerable scatter 
among the a s (M^) values from all the observables. In each case the average over 
the a s (M^) values, and corresponding r.m.s. deviation, are shown in Table 1. The 
r.m.s. deviations are comparable with those resulting from choice of the physical and 
experimentally-optimised scales, implying that higher-order effects contribute roughly 
equally in all of these procedures. 

An approach for estimating higher-order perturbative contributions to, as well as 
the sum of, perturbative QCD series is based on Pade Approximants (PA). The PA 
[N/M] to the series: 

S = So + S 1X + S 2 x 2 + ... + S N+M x N+M (2) 

is defined [14] 



where N and M are integers such that N > and M > 0, and 

[N/M] = S + 0(x N+M+1 ). (4) 

The coefficients a n (0 < i < N) and bj (I < j < M) are obtained by multiplying 
Eq. 4 by the denominator of Eq. 3 and equating coefficients of like powers of x. By 
consideration of the terms of (^+ M+1 )one can obtain an estimate of the coefficient 
iSjv+m+i- Furthermore, for an asymptotic series [N/M] can be taken to be an estimate 
of the sum (PS) of the series to all orders. 
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[N/M] 



ao + a\x + a>2X + . . . + aj^x 
1 + b\x + b 2 x 2 + ... + b M x M 



(3) 



In the case of hadronic event shape observables the PA [0/1] can be defined for 
the series Eq. (1) and, for each bin of each observable, was used to derive an estimate 
of the coefficient C of the O(a^) term [15]. The PA prediction for C was added to 
the exact 0(a 2 ) calculation to obtain an estimate of the series to O(a^). For each 
observable the calculation was fitted to the data [5] using // = Mz, and the resulting 
(Aff ) values are shown in Fig. 2(a). The O(a^) estimate does not provide a good 
fit to the Bt data [15] and this observable is excluded from further discussion. In 
each case the a s (Af|) value derived using the O(a^) estimate is lower than that 
derived using the 0(a 2 ) calculation, which is expected since C is positive [15], and the 
O(af) a s (M^) value lies near the lower bound given by the scale uncertainty on the 
0(a 2 ) result. To the extent that the PA O(a^) estimate is accurate, this implies that 
the renormalisation scale uncertainty assigned to the 0(a 2 ) a s (M^) value from each 
observable is a reasonable estimate of the effect of the missing O(af) contribution. 
The average and r.m.s. deviation of the 14 a s (Af|) values are listed in Table 1. 
The scatter is noticeably smaller than in any of the 0(a 2 ) cases, implying that the 
Pade method provides at least a partial approximation of higher-order perturbative 
QCD contributions to event shape observables. 

Finally, the PS [0/1] was used as an estimate of the sum of the asymptotic series 
and a s (Af|) was extracted by comparison with the data in a similar manner [15]. 
The a s (Af|) values are shown in Fig. 2(b). Typically, for each observable, the PS 
a s (Af|) value is close to the PA O(af) value. Again the fit to Bx is very poor [15]. 
The average and r.m.s. deviation over the set of 14 a s (Af|) values are listed in Table 
1. Though the average value is close to that obtained using the PA O(a^) estimate, 
the r.m.s. deviation is somewhat larger, implying that the PS [0/1] provides a poorer 
estimate of the sum of the series than the PA [0/1] estimate to O(a^). 
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In summary, a s (Af|) has been determined by fitting O(o^) QCD predictions of 
f5 hadronic event shape observables to e + e - annihilation data at the Z° resonance 
collected by the SLD experiment. Five prescriptions for choosing the renormalisation 
scale were used, namely the physical, experimentally-optimised, PMS-, FAC- and 
BLM-optimised scales. Though the average a s (Af|) value, taken over all the 
observables, differs among these five procedures, the scatter among the a s (Af|) values 
from different observables is equally large in each case, the r.m.s. deviation being 
about 0.008, implying that these specific renormalisation scale choices do not offer any 
numerical advantage in terms of the accuracy of O(o^) perturbative QCD predictions 
of e + e - event shapes. 

If Pade Approximants are used to estimate the O(al) terms the scatter among the 
a s (Af|) values from different observables is reduced to ±0.0035. This is comparable 
with the combined experimental error and hadronisation uncertainty on a single 
observable measured at Q = Mz [5]. Since the accuracy of the Pade Approximant 
method can only be verified a posteriori, exact calculation of the O(af) terms in order 
to confirm these results is extremely desirable. 
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Procedure 


a s (Ml) 


Physical scale 


0.1265 ± 0.0076 


Exp. opt. scale 


1173 ± 0071 


PMS scale 


1123 ± 0079 


FAC scale 


0.1123 ± 0.0080 


BLM scale 


0.1088 ± 0.0075 


Pade 0(a 3 s ) 


0.1147 ± 0.0035 


Pade sum 


0.1148 ± 0.0052 



Table 1: Mean and r.m.s. a s (Af|) values determined using different theoretical 
procedures. 

Figure Captions 

FIG. 1. Values of a s (Af|) from QCD fits to the data using: (a) physical (solid 
circles), and experimentally-optimised (open circles) scales; (b) PMS- (solid circles), 
FAC- (solid triangles), and BLM- (open squares) optimised scales. In all cases only 
experimental error bars are shown. For each observable the shaded region indicates 
the total uncertainty estimated in Ref. [5], dominated by the contribution from wide 
variation of the renormalisation scale. 

FIG. 2. Values of a s (Af|) from QCD fits to the data using: (a) PA 0(0^) estimate 
(squares); (b) Pade sum (PS) (crosses). In all cases only experimental error bars 
are shown. For each observable the shaded region indicates the total uncertainty 
estimated in Ref. [5], dominated by the contribution from wide variation of the 
renormalisation scale. 
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